Master equation approach for interacting slow- and stationary-light polaritons 
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A master equation approach for the description of dark-state polaritons in coherently driven 
atomic media is presented. This technique provides a description of light-matter interactions under 
conditions of electromagnetically induced transparency (EIT) that is well suited for the treatment 
of polariton losses. The master equation approach allows us to describe general polariton-polariton 
interactions that may be conservative, dissipative or a mixture of both. In particular, it enables us 
to study dissipation-induced correlations as a means for the creation of strongly correlated polariton 
systems. Our technique reveals a loss mechanism for stationary-light polaritons that has not been 
discussed so far. We find that polariton losses in level configurations with non-degenerate ground 
states can be a multiple of those in level schemes with degenerate ground states. 

PACS numbers: 42.50.Gy,32.80.Qk,42.50.Xa,42.65.-k 



O 



Oh' 
I 



(N 
> 

in 

00 

^, 

o 
o 



% 



INTRODUCTION 



Photons are ideal carriers for quantum information 
over long distances. This is due to the large propaga- 
tion speed of light and the fact that photons in free 
space do not interact with each other. On the other 
hand, the generation of highly entangled light fields and 
the realization of photon gates requires strong photon- 
photon interactions [l| . Nonlinear media can mediate an 
effective interaction between photons, but the strength 
of this induced coupling is usually weak. Thus the re- 
alization of strong photon-photon interactions is a ma- 
jor challenge in quantum information science. Similarly, 
strong photon-photon interactions are a key requirement 
for quantum optical irnplementations of highly correlated 
many-body systems [2|. A substantial research effort 
|3l-[l4l| is currently devoted to these systems where com- 
bined excitations of light and matter, i.e. polaritons, 
reproduce the dynamics of bosons with tunable mass 
and different interaction types. Several effects in corre- 
lated many-body systems were considered, including the 
realization of Bose- Hubbard models 0t3 [3 ^ quantum 
transport (8|j_9[, nonlinear effects in driven dissipative 
systems [3, [l5 1 , Bose-Einstein condensatio n fill , a nd the 
realization of a Tonks- Girardeau gas [10, [l2|, [Tj] of po- 
laritons. 



Dark-state polaritons |16l - ll8l ] represent bosonic quasi- 
particles that arise in light-matter interactions under 
conditions of electromagnetically induced transparency 
(EIT) [1^. The generic EIT scheme consists of a gas 
comprised of three-level atoms in A configuration that are 
driven by a strong control field and a weak probe field on 
separate transitions. EIT gives rise to a multitude of in- 



triguing effects like the slowing and stopping of li ght 20 - 
22 1 , the coherent storage and retrieval of light [la, [13, [23 - 
25 1 and stationary light [ill. [26l - l33l |. Of particular rel- 
evance for quantum optical realizations of strongly cor- 
related many-body systems is stationary light where the 
polaritons experience a quadratic dispersion relation like 
bosons in free space. Moreover, EIT in atomic four-level 
systems gives rise to a strongly enhanced Kerr nonlin- 



earity for the probe fields [3J-|40|. This effect results in 
a two-particle contact interaction between the polaritons 
that can be conservative [l^, |4l| or dissipative [13 • 

Experimental setups that show great promise for the 
realization of strongly correlated polariton systems are 
arrays of coupled microcavities doped with emitters [42|- 
l46j or optical fibers that couple to atoms [43, 13 • I'^ 
general, strong light-matter interactions require the con- 
finement of light to small interaction volumes. Here we 
consider the experimental setup [47| shown in Fig. |TJ 
where photons and atoms are simultaneously confined 
to the hollow core of a photonic-crystal fiber, and the 
level scheme of each atom is shown in Fig. |2j Since the 
light-guiding core of the optical fiber is of the same order 
of magnitude as the optical wavelength, the fiber repre- 
sents a one-dimensional waveguide for the optical fields. 
Note that a second potential realization is comprised of 
the experimental setup in |48| , where multi-color evanes- 
cent light fields surrounding an optical nanofiber couple 
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FIG. 1. (Color online) Considered setup of A*' atoms confined 
to an interaction volume of length L and transverse area A. 
f2± are the Rabi frequencies of the classical control fields, and 
E± are the quantum probe fields. A typical intensity profile 
of the quantum fields is shown for the case of stationary light 
where the control fields have the same intensity {Q.+ = $!_). 
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of MaxwcU-Bloch equations and find excellent agreement. 
The full master equation including the most general form 
of polariton-polariton interactions is covered in Sec. lIVBl 
and the mapping to the dissipative Lieb-Liniger model is 
outlined in Sec. IIVCI Finally, in Sec. IIVDI we compare 
the advantages and disadvantages between the atomic 
level schemes in Figs. [2] and [5] that both give rise to the 
same master equation for dark-state polaritons. 



FIG. 2. (Color online) Atomic level scheme, ytj is the full 
decay rate on the |i) <->■ |j) transition, 5 and A label the de- 
tuning of the probe fields with states |3) and |4), respectively, 
and e is the two-photon detuning. 



to atoms trapped in an optical lattice. 

Existing descriptions [la, ll3 of dark-state polaritons 
in EIT systems are based on a Heisenberg-Langevin ap- 
proach for the polariton field operator. Here we present a 
different approach and derive a master equation for the 
reduced density operator of dark-state polaritons. The 
master equation technique facilitates the treatment of 
polariton losses and allows one to account for general 
polariton-polariton interactions that may be conserva- 
tive, dissipative or a mixture of both. This is an im- 
portant achievement since it opens up the possibility to 
study dissipation-induced correlations j49l ISOj in polari- 
ton systems. Dissipative polariton-polariton interactions 
are promising in the quest for highly correlated systems 
since they can be considerably stronger than their conser- 
vative counterparts [l3|. Second, our method reveals an 
additional loss term for stationary-light polaritons whose 
importance depends on the structure of the atomic level 
scheme and that was not discussed in the literature yet. 

This paper is organized as follows. In Sec. |TT] we set 
up a master equation for the atoms interacting with the 
quantized probe and classical control fields inside the ID 
waveguide, see Figs. [1] and [2l We then transfer the orig- 
inal master equation into a master equation solely for 
dark-state polaritons. This process is detailed in Sec. IIIII 
and consists of four steps. First, we show that the en- 
tire problem can be described in terms of bosonic quasi- 
particles if the number of atoms is much larger than the 
number of probe field photons, see Sec. IIII Al The con- 
cept of dark-state polaritons is introduced in Sec. IIIIBi 
and the formulation of the original master equation of 
Sec.|TT]in terms of dark-state polaritons and all other ex- 
citations is presented in Sec. IIII CI In the final step of 
the derivation we trace out all excitations except for the 
dark-state polaritons, see Sec. IIIIDI The master equa- 
tion for dark-state polaritons under conditions of station- 
ary light and for arbitrary (conservative or dissipative) 
polariton-polariton interactions is presented in Sec. IIVI 
Here we summarize all conditions that grant the valid- 
ity of our approach. The special case where the level 
scheme in Fig. [5] is reduced to the A subsystem is dis- 
cussed in Sec. IIV Al We compare the predictions of our 
master equation to the results of a numerical integration 



II. DESCRIPTION OF THE SYSTEM 

We start with a more detailed description of our one- 
dimensional model shown in Figs. [T] and [21 Each of the N 
atoms interacts with control and probe fields denoted by 
il± and E±, respectively. The control fields of frequency 
ujc are treated classically and 51+ (fl-) labels the Rabi 
frequency of the control field propagating in the posi- 
tive (negative) z direction. In addition, we assume that 
the control fields are spatially homogeneous and that the 
Rabi frequencies Q± are real. The probe fields E+ and 
E- are quantum fields that propagate in the positive and 
negative z direction, respectively. They are defined as 



E±{z) = ^a±k^+ke 



i{±k^ + k)z 



(1) 



where a±k^+k are photon annihilation operators of a 
mode with frequency uj±k^+k and kc (— fcc) is the wave 
number of the control field ft^ {^-)- We assume that the 
wave numbers k satisfy |fc| ^ kc which implies that the 
envelope of the quantum fields varies slowly on a length- 
scale defined by the wavelength of the optical fields. 
We model the time evolution of the atoms and the 



quantized probe fields by a master equation 51| for their 
density operator g, 



g = --[H,g] 



C^g, 



(2) 



where the system Hamiltonian H = Hq + H\ + i^NL 
is comprised of three parts. Hq describes the free time 
evolution of the atoms and the probe fields, and H\ ac- 
counts for the interaction of the probe and control fields 
with the A-subsystem formed by states |1), |2) and |3). 
On the other hand, i?NL accounts for the coupling of the 
probe fields to the |2) <-> |4) transition and results in a 
nonlinear coupling between probe field photons [3J-|40|. 
In a rotating frame that removes the time-dependence of 
the classical laser fields, Hq, Ha and i^NL are given by 



Ho = 



- hJ2i^P - ^k^+k)al^^^ak,+k (3) 
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H^L = - hg2 J2 Si2 E+ (z^) + E^ (z^) 



h.c. . 



(5) 



^L=l 



Af^) 



Here A''^ and S^'^ are projection and transition operators 
of atom fi at position z^, respectively, 



^^ = lv)(vL 



s. 



M 



Iv)OmI (*^j-), (6) 



ujp is the central frequency of the probe pulse, and gi 
((72) is the single-photon Rabi frequency on the |3) <-^ |1) 
(1 4) o 1 2)) transition. The detuning of the probe field 
with respect to the transition |3) O- |1) (|4) o |2)) is 
labeled by S (A), and e is the two-photon detuning, 



S^lVp-UJsi, A^UJp-UJ4,2, 



Up -Uc- UJ2. (7) 



Here the energy of level \i) is fiuii (we set oji ~ 0) and 
transition frequencies are denoted by w^ = oj^ — ojj . The 
term C-yg in Eq. ^ accounts for spontaneous emission 
from states |3) and |4), 
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where 7^^ is the full decay rate on the transition \i) <-)• \j) 
(see Fig.[T]). Finally, we introduce the parameter 



Aw 



LJ2 



(9) 



are in state |1) and that the total number of photons is 
much smaller than the number of atoms TV. This assump- 
tion allows us to study the system dynamics entirely in 
terms of independent bosonic quasi-particle excitations, 
see Sec. lIH Al A second simplification is made possible by 
the concept of dark-state polaritons [1^, [ill introduced 
in Sec. IHIBI Dark-state polaritons are bosonic quasi- 
particlcs that decay only indirectly via the coupling to 
other bosonic modes that are termed bath excitations. 
In the so-called slow-light regime, this coupling is much 
slower than the decay of the bath excitations which is 
of the order of the decay rates of the excited states |3) 
and 1 4). The existence of these two different time scales 
enables us to derive a Markovian master equation for the 
reduced density operator of the dark-state polaritons as 
outlined in Sec. IHIDl Throughout this Section, all tech- 
nical details and lengthy definitions are moved to the 
Appendix. 



A. BOSONIZATION 

Here we show that the system described in SecHIlcan 
be mapped to a much simpler system if almost all atoms 
are in state 1 1) and if the total number of photons is much 
smaller than the number of atoms A^. We begin with the 
description of a simple system that consists of M bosonic 
modes. The annihilation operator of mode i is given by 
Oi {ie {1,..., M}) and O = {d, . . . , Om} denotes the 
set of all operators that obey the commutation relations 



4,oI 



Hi J), 



o^,o. 



(10) 



If |0) is the vacuum state of the system, it follows that 



M 

|{ni,...,nM}>=n^=f(0|)"lO) (11) 



which describes the frequency difference between the 
probe and control fields. Note that Aw is practically 
equal to the frequency splitting W2 between the ground 
states |1) and |2) if the two-photon detuning e is small. 
Next we outline the approach we developed to reduce the 
master equation ^ for the atoms and quantized probe 
fields into a master equation solely for dark-state polari- 
tons [la] , formed by collective excitations of photons and 
atoms. 



III. MASTER EQUATION FOR DARK-STATE 
POLARITONS: DERIVATION 



Here we show that the master equation ([2]) can be sim- 
plified considerably if we assume that almost all atoms 



is a normalized state with rii excitations in mode i and 
a total number of X)i=i '^i excitations. Furthermore, we 
note that the total state space of M bosonic modes is the 
tensor product of the state spaces Hi associated with the 
individual modes, 



Hose =Hi(E)H2<S)...<S)Hm- 



(12) 



Next we show how the system of Sec.|TT]can be mapped 
to this simple model outlined in Eqs. ([TU1) - P^ . First we 
define a vacuum state where all probe field modes are 
empty and all atoms are in state |1), 



|0) = |{0}phot;li,...,iiv>. 



(13) 



Second, we define the following operators {m € Z) 

Ak = ak^+k sin Lp + a^k,+k cos <^, 
Dk = flfee+fc COS ip - a^k,+k sin (/?, 

AT 



^w^^Es! 



(/j) -i(nifce+fe)Zf, 



TV 



u=l 



•3 ^ ' 



him) 



1 V^ c'(m) 



(m) -i(rnfce+fc)z^ 



(14) 
(15) 

(16) 
(17) 
(18) 



/.= ! 



where k = n2n/L [n e Z) and A^ (£>fc) is a sum (dif- 
ference) of two counter-propagating probe field modes. 
Since a±k^+k are photon annihilation operators, it fol- 
lows that Ak and Dk obey bosonic commutation rela- 
tions. The angle ip depends on the relative strength of 
the Rabi frequencies ^+ and Vt- and is defined by |28l.l29j 



sin(^ 



r^i 



nl + n''_ 



COS(p = 



n^ 



nl 



91 



(19) 



The operator Xk (m) describes a collective spin coherence 
that is slowly oscillating for to = and fast oscillating 
for m 7^ 0. The operators Hl{m) and lli'm) create an 
excitation in the excited states |3) and |4), respectively. 
Next we show that the operators defined in Eqs. (fro|) - (fT5)) 
obey bosonic commutation relations for all wave numbers 
k and all to € Z if almost all atoms are in state |1). As 
an example, we discuss the commutation relations for 
Xk{m). Within a manifold with fixed m, we find 



1 '^ 

[X,(TO),Xt(TO)]=^5](<-4^^ 



N 



i{m) _ M\ ^i{p-k)z 



-E 



J{p-k)z^ 



(20) 



(21) 



L 



]- f dze'^P'^^'' = d{k,p), 



where we set A\1' sa 1 and A22' ~ since almost all 



1^^^ « 1 and A^f^^ 
atoms are in state |1). Furthermore, we employed that 
the mean distance between atoms is much smaller than 
l/|fc| for all relevant wavenumbers k contributing to the 
slowly varying envelopes of the control fields. In the case 
m ^ n, wc find 



1 ^ 
[Xk{m),Xl{n)]^-J2 



^i[{p-k) + (n~m)kc]Zf^ 



(22) 



In contrast to Eq. (j20[) , the sum cannot be converted into 
an integral since the mean spacing between the atoms is 
much larger than 1/fcc for realistic densities, where kc is 



the wavenumber of an optical transition. However, the 
sum in Eq. (|22p represents the average of random num- 
bers on the unit circle in the complex plain, which is zero 
in the limit N ^ 00. We can thus set [Xk{m), X^{n)] « 
for m ^ n. Since [Xk{m),Xp{n)] = 0, it follows that the 
operators Xk{m) obey bosonic commutation relations if 
almost all atoms are in state |1), and corrections scale 
with 1/A^. The same result is found for Hk{m) and 
Ik{m). Furthermore, we point out that Xk{m), Hk{m) 
and Ik{iTi) describe independent excitations up to cor- 
rections that scale with l/iV, i.e., [Xk{m),H^{n)] w 0, 
[Xfe(TO),/pt(n)] « and [Jfe(TO),iJt(n)] « 0. 

In summary, we can introduce the set of bosonic oper- 
ators 

O ^ {Ak,Dk,Xk{m),Hkim.),Ikim)}. (23) 

if almost all atoms are in state |1). This condition can be 
met if the total state space "Htot of the system is restricted 
to the subspace 



■Hfe = Span< |{ni, . . . ,71^}) | ^n.j < iV ^ 



(24) 



that is spanned by states with much less excitations than 
number of atoms A^. From a physical point of view, the 
system dynamics will be restricted to this subspace if the 
number of photons is much smaller than the number of 
atoms, and if initially almost all atoms are in state |1). 
In Appendix \^ we show that the Hamiltonian H and 
the decay term £^g can be expressed entirely in terms 
of the bosonic operators in Eq. 1^^ if the state space 
is restricted to ?^fe- We denote the density operator in 
^FE by p, and the master equation ^ in 'Hfe can be 
written as 






C^Q, 



(25) 



where H = Hq + Ha + -^nl- The operators Hq, H\, 
Hnl and £^g arc defined in Eqs. (|A12|) . (jA14|) . (|A15|) 
and (jA16|) . respectively. These operators approximate 
their counterparts without tilde in the subspace "Hfe and 
are comprised of the bosonic operators in Eq. 



B. DARK-STATE POLARITONS 

The interaction Hamiltonian H\ of the A-subsystem 
has the important property that a certain class of its 
eigenstates are so-called dark states \D). These states 
are called dark since they do not contain a contribution 
of the excited state |3) and are thus immune against spon- 
taneous emission. A simple example for a dark state is 
given by 



i^>=^'iio), 



(26) 



where the unique definition of the operator tpk is J29| 

tpk = AkC0s9 — XkSinO. (27) 



In this equation, Ak is a superposition of two counter- 
propagating probe field modes, and Xk = ^fc(O) de- 
scribes a slowly varying collective spin coherence. The 
mixing angle 9 determines the weight of the photonic 
(Ak) and atomic (Xfe) components contributing to ipk 
and is defined as 0, M, IM M HI] 



m, „_, \/"i + » 


(28) 


JIq ^0 


'Ngj + nl + ni. 


(29) 



nn 



A short calculation shows that \D) in Eq. (j26|) is an 
eigenstate of Ha in Eq. ^ with eigenvalue zero, i.e., 
Ha\D) =^ 0. Furthermore, Eq. (gT)) implies that \D) does 
not contain a contribution of the excited state |3). It 
follows that state \D) is indeed a dark state of the in- 
teraction Hamiltonian Ha- Note, however, that \D) is 
not an eigenstate of the remaining parts Hq and i?NL 
of the full Hamiltonian, and these terms give rise to a 
non-trivial time evolution of the dark-state polaritons. 

The results of Sec. IIII Al and Eq. (gT]) imply that the op- 
erators ipk obey bosonic commutation relations in Hfe, 



[V'fc,V'p] ~^{k,p), 



(30) 



and the quasi-particles associated with these excitations 
are termed dark-state polaritons. It follows from Eq. ((50)) 
that 



1^) 




(31) 



represents a normalized dark state with ^^. n^ excita- 
tions. We emphasize that the photonic part A^ of ipk 
contains a pair of counter-propagating probe field modes 
that are grouped around the wavenumbers of the con- 
trol fields ±fcc rather than the mean wavenumbers ±fcp 
of the probe fields, see Eq. ((Ti)) . Note that the states in 
Eq. (pri) would not be true dark states if the probe field 
modes in A^ were grouped around ifcp and if both i7+ 
and 51- were different from zero. 



BOSONIZATION WITH DARK-STATE 
POLARITONS 



We have shown in Sec. IIII Al that the master equa- 
tion ([2]) can be formulated in terms of bosonic modes if 
the system dynamics is restricted to the subspace Hfe- 
Here we restate this model in terms of long-lived dark- 
state polaritons introduced in Sec. IIIIBI With the defi- 
nition of the bright-state polariton [11[ 



Ak sin 9 ~\- Xk cos 9 



(32) 



and the inverse relations of Eqs. (P?]) and ([5^. we replace 
the operators Ak and Xk{0) in O [see Eq. ((23)) ] by ^k and 
(pk- This allows us to write the subspace "Hfe in Eq- (EH) 



as the tensor product of the state space Hs of dark-state 
polaritons and the state space "Hb of all other modes 
termed bath excitations, 



■Hfe ='Hs<E}Hb- 



(33) 



The partition of all bosonic modes into dark-state po- 
laritons and bath excitations is motivated by our aim 
to derive a master equation for the long-lived dark-state 
polaritons only, see Sec. IIIIDI 

In the following, we assume that the Rabi frequencies 
of the control fields are identical and set 



ilr — J M 



fi_. 



(34) 



With this choice. Ha gives rise to the stationary light 
phenomenon [2a, l23| that allows us to trap the probe 
field inside the medium. Note that any other choice of the 
Rabi frequencies il,± can be treated within the formalism 
introduced here, and in these cases the calculation follows 
exactly the same route as detailed below. If the operators 
Ak and Xk{0) in the master equation (^5)) are replaced 
by Ipk and pk, we obtain (see Appendix IBJ) 

g = -^[Hs,g]-^[V,g] + CBg, (35) 

where 

Hs = -h{sm 9e + cos^ 9) ^ VlV'fc (36) 

k 

describes the free time evolution of the dark-state polari- 
tons. Here we choose a small two-photon detuning 



e = — cot" 9/S.lli 



(37) 



such that Ha, = 0. The dynamics of the bath excitations 
is governed by the Liouvillian 



C-eg 



\H 



B,g\ 



/:f 



g, 



(38) 



where iJe accounts for the unitary time evolution of the 
bath modes and the decay of bath excitations is deter- 
mined by C\ g. The interaction between dark-state po- 
laritons and other excitations in 'Hb is described by the 
interaction Hamiltonian V. The definitions oi Hb, C\ q 
and V are provided in Appendix [BJ 

(B\ 

The decay term C\ g results in a finite lifetime of 
the bath excitations that is of the order of the lifetimes 
of the excited states |3) and |4). On the other hand, 
the dark-state polaritons decay only indirectly via the 
coupling to bath excitations mediated by V. In the slow 
light limit, this coupling is much slower than the decay 
of the bath excitations. This existence of two different 
time scales opens up the possibility to derive a Markovian 
master equation for the dark-state polaritons alone, and 
this procedure is outlined in the next Section UlIDI 



D. ELIMINATION OF THE BATH 



where the rising and lowering parts of V are defined as 



In the previous Sections IIII AfllTC] we achieved to 
transform the initial master equation ^ within the sub- 
space Hfe into a master equation for long-lived dark- 
state polaritons and fast-decaying bath excitations. Here 
we are especially interested in the quantum state gu of 
the dark-state polaritons that is obtained from g by a 
partial trace over all excitations except for the dark state 
polaritons. We derive the corresponding master equa- 
tion for pp from Eq. psp via projection operator tech- 
niques |5l| and assume that the initial state of the system 
factorizes into a product of the initial polariton state qd 
and the vacuum state qb of the bath modes, 



g{t = 0) = £iD «> Qb- 



(39) 



Furthermore, we employ the Born-Markov approxima- 
tion [sii and obtain 



gh = -S{gu) + h.c. 



(40) 



where 



Sign) 



1 



cItTtb { [V, e^^^Vgj, ® g^]] ■ (41) 



l/(+)=^Bt5,, V^-^^Y.BA 



(43) 



respectively. In this equation, Si represents a system 
operator comprised of dark-state polaritons, and Bi is a 
bath operator. Since we assume that the bath is initially 
in its vacuum state, we have V^'g-B = which allows 
us to replace V gB by l/^+^gie in Eq. (|¥T|) . In addition, 
the second interaction Hamiltonian V in Eq. ([^T|) can 
be replaced by V^~^ since the contribution of T/(+) is 
negligible, see Appendix [C] We thus arrive at 



5(gD) = ^ I drTvB { [F(-),e^-^V^W 



gD <x) gB 



OO 



(44) 



»J 



X [Sls^gxi - SjgxiS'^^ , 



and the evaluation of the bath correlation functions 
Trel-Bi exp[£B''"]SJ£'B} is presented in Appendix [U] The 
final result for the master equation (PO)) is discussed in 
the next Section HVl 



The application of the Born-Markov approximation re- 
quires that the coupling of the dark-state polaritons to 
bath excitations is sufficiently weak and in particular 
small as compared to the decay rate of bath excitations. 
Conditions for the validity of the Born-Markov approx- 
imation as well as the assumption in Eq. p9|) are dis- 
cussed in Sec. IIVI and Appendix O In order to outline 
the evaluation of Eq. (|4T]) , we write V as 



(42) 



F = y(+) + y(-), 



IV. MASTER EQUATION FOR DARK-STATE 
POLARITONS: RESULTS 

The master equation ^ in Sec. [Tl] describes the in- 
teraction of classical control and quantized probe fields 
with A^ atoms. In Sec. IIIII we demonstrated that this 
master equation can be converted into a master equation 
for the reduced density operator g^ of dark-state polari- 
tons. For equally strong control fields [see Eq. (j34p ] and 
in the slow-light limit cos'^ 6' <C 1 [see Eq. ((29)) ]. we obtain 



PD = - T[Hnd, Pd] - T [i?3, Qn] - T [i?4, gn] - 7^ COS^ 6* ^ CJ^ (V'iV'fcfi'D + PDV'lV'fc " 2i'ugD^l 



(45) 



2 „„„2 



2^2 



Aw cos 



.gj 742/2 

A^+7|2/4 



^y^ yPli^kg-D + goi'lipk - 2^pkgDipl 

k 

cos^ ^ yiPl'>Pl_ptpq'>Pk~qgD + gi)i^li^l_p'>Pqipk-q - 2-4>q->pk-qgi)tPlipi_pj , 



k,p,q 



where 



i?nd = -h—^Auj^ cos^ 9j2'^i'^k, 



(46) 



^3 = -h-^ COS^ J2 '^fe^i^fc' 



Hi 



hAggl 



■7I2/4 



COS^ e Y^ ^l^i_p^q^, 



k—q- 



(47) 
(48) 



k,p,q 



Here Aw 



LOc is the frequency difference between 



the probe and control fields, 



Afi 



A - cot^ 6* Aw, 



(49) 



and r = 731 + 732 is the full decay rate of state |3). 

Next we summarize the conditions under which 
Eq. (|^5|) holds. First of all, we note that the two-photon 
detuning e is constrained by Eq. ((37)) , and it was assumed 
that fio defined in Eq. ^^ is large as compared to the 
decay rates of the excited states and the detuning with 
state 1 3), 



^o»r,742,|5|. 



(50) 



A key assumption in the derivation of Eq. (P5| is that 
the number of atoms is much larger than the number of 
photons, and that initially almost all atoms are in state 
|1). This condition is a prerequisite for the bosonization 
described in Sec. IIII Al Furthermore, it was assumed in 
Sec. IIIIDJ that initially all bath modes are in the vacuum 
state, and that the initial density operator of the total 
system factorizes, see Eq. ([39)) . These conditions can be 
met if the initial dark-state polariton state is prepared 
via the slowing and stopping of a probe pulse [16, 17, 23- 
m, m US mi- in tns case, the initial state of the dark- 
state polaritons is a slowly-varying spin coherence, and 
all other modes are in the vacuum state. The regime 
of stationary light can then be entered if the counter- 
propagating control fields arc adiabatically switched on. 
The derivation of Eq. ps)) relies on the validity of the 
Born-Markov approximation. The Born approximation 
employed in Sec. IIII DJ reauires that the coupling between 
the dark-state polaritons and other excitations is suffi- 
ciently weak and holds in the slow-light regime where 
cos^ 9 <C 1. On the other hand, the Markov approxi- 
mation relies on the existence of two very different time 
scales. The fast time scale is represented by the bath 
correlation times T-q and is of the order of the lifetimes 
of the excited states |3) and |4). On the other hand, the 
slow time scale Tg is given by the typical evolution time 
of dark-state polaritons. The condition Ts >• Tb which 
justifies the Markov approximation is fulfilled provided 
that the following inequalities hold. 



cos^ Oc^k'^. 

"0 
IQgl cos2 6liVph 

7I2 



21,2 



<1, 



<1. 



2\5\ cos^ Oc^k^ 
W 



<1, 



(51) 
(52) 



Here c is the speed of light in the fiber, A^ph is the number 
of photons in the pulse and fcmax > \k\ is the maximum 
of all occupied wave numbers. In addition, we empha- 
size that the Markov approximation is only possible if 
we assume that the fast ground-state coherences Xk{m) 
for m ^ [see Eq. p^ ] are washed out due to the atomic 
motion on a timescale comparable to the lifetime of the 
excited state |3), see Appendix [Cj Note that this condi- 
tion is not required in the case of the level scheme dis- 
cussed in Sec. IIVDI If Fpo denotes the decay rate of the 



fast oscillating spin coherences, the validity of the master 
equation (|45p requires 



(Oo cos 61)2 
Epo 



<r. 



(53) 



If this condition is not fulfilled but if Fpo is of the order 
of the decay rates of the excited states, then the gen- 
eral structure of Eq. (|45[) remains the same, but the pre- 
factors of iJ3, i^nd and the decay terms proportional to 
F will be different. 

Finally, we note that the Hamiltonian i7nd in Eq. pS)) 
represents a constant energy shift of the polariton exci- 
tations. This term is only present if Aw 7^ and hence if 
the ground states |1) and |2) are non-degenerate. In the 
following, we will work in an interaction picture with re- 
spect to iJjid , and the representation of the master equa- 
tion in this rotating frame can be obtained if i^nd is omit- 
ted in Eq. (|i5)) . 



A. STATIONARY LIGHT 

In this Section we focus on the phenomenon of sta- 
tionary light [23, [23[ that arises from the interaction of 
the probe and control fields with the A system formed 
by states |1), |2) and |3). Formally, the reduction of the 
master equation (|45)) to this case is accomplished if the 
coupling constant 52 is set equal to zero. In the follow- 
ing, we formulate the master equation in terms of the 
operator 



^w-:^E 



^/L■ 



ikz I 



(54) 



It follows from Eq. (j30p that -i/! is a bosonic field operator 
that obeys the commutation relations 

mz),i:Hz')]=S{z^z'), [V'(z),V(^')]=0. (55) 

The master equation (|15|) for (72 = reads 



BB 



t\-^3, go] + CigB + C2Bt>, 



(56) 



where H^ is defined in Eq. (|T7)) and shows that the po- 
laritons experience a quadratic dispersion relation. With 
the deffiiition ()54[). this Hamiltonian can be written in 
the form of a kinetic energy term. 



i?.^ = 



2to 



eff 



dzd:,'ip''dzip, 



where 



TOcff 



hnl 



2Sc^ cos2 9 
is the effective mass of the polaritons. 



(57) 



(58) 
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FIG. 3. (Color online) Decay of the fc = polariton mode 
for non-degenerate ground states |1) and |2). The solid line 
shows the decay according to Eq. (|62ll . the corresponding re- 
sult from a numerical integration of Maxwell-Bloch equations 
is indicated by red dots. The parameters are Auj = 150r, 
cos^6l = 2.5 X 10'", ^0 = 90r, (5 = and £ is fixed by 
Eq. (|37|l . The initial spin coherence R21 (see Appendix ID)) 
has a Gaussian shape -R21 oc exp[— (z — zo)'^ /{2a'^)] with 
a = 24.3rc/f2o) s-nd the control fields are ramped up accord- 
ing to ndt) = 0.5[1 -I- tanh(f - lO/F)] x T. The initial value 
(V'oV'o)init refers to {tpQipo) at i = 20/r where the control fields 
have attained their maximal value of Qc ~ E. 



The terms Cigu and £2£'d in Eq. ([55)) describe polari- 
ton losses and are defined as 



Cigu = ~^^c^ cos' OV[d,^P] 



^■29-0 = 
respectively, where 

L 



2ni 



Aw2cos2 6'I?[V'], 



(59) 
(60) 



V[X] = I Az{X^XgY> + QYiX'^X - 2X^0^^) (61) 



is a dissipator in Lindblad form J5lJ for an operator X . 
The term £i£id arises due to the coupling between dark- 
state polaritons and the difference mode Dk- Since the 
decay rate of the individual modes increases quadrati- 
cally with the wave number, CiQd does not result in an 
exponential damping but leads to diffusion [23, [2^ • On 
the contrary, L2QT) leads to identical decay rates for all 
modes. This term stems from the coupling between dark- 
state polaritons ^k and bright polaritons (jy^. Since this 
loss mechanism has not been discussed in the literature 
yet, we investigate it in more detail here. First of all, 
we note that C2Qt> is proportional to Aw^, and Aw prac- 
tically coincides with the splitting of the ground states 
|1) and 1 2) [see Eq. ^]. An important consequence of 
C2Qt> is that in contrast to EIT, dark-state polaritons in 
the fc = mode decay under the conditions of stationary 
light provided that the ground states |1) and |2) are non- 
dcgcncratc. The decay of the mean number of dark-state 



FIG. 4. (Color online) The number of dark-state polaritons 
is shown as a function of time. The results of Eqs. (|64|) and 
Eq. (|65l) are represented by the solid and the dotted line, re- 
spectively. The good agreement indicates that the polariton 
losses are described correctly by the master equation (|56|l . 
The number of polaritons follows the dashed line if the addi- 
tional decay term £2^10 in Eq. (|56|) is omitted. The parame- 
ters are the same as in Fig. [31 



polaritons in the A: = mode can be calculated from 
Eq. (|56p and is given by 

dtii^l^o) = -T^Ac^2^os2 0(^.J^o). (62) 

"0 

The accuracy of this result can be confirmed numerically 
if (V'oV'o) is evaluated via MaxwcU-Bloch equations for 
classical probe fields (see Appendix |D]). The result is 
shown in Fig. [3l where the solid line corresponds to the 
exponential decay according to Eq. (|62p . The dotted line 
represents (V'qV'o) obtained from the numerical integra- 
tion of MaxwcU-Bloch equations and is in perfect agree- 
ment with the predictions of the master equation (|56p . 

Next we compare the impact of the loss terms Ciq-q 
and L2Qt>- Equations ([59)1 and ([5(I[) imply that £2£'d 
becomes comparable to £igD if |Aw| is of the order of 
cfcmax • Since cfcmax ~ | Aw | -|- CGk whcrc (Jk is the width 
of the polariton pulse in k space, the two loss terms are 
comparable if 



I Aw I > cuk- 



(63) 



On the other hand, the width c^ can be estimated to be 
of the order of 2tt/L, where L is the length of the system. 
The inequality (j63p thus implies that the impact of C2QD 
is of the same order of CiQd if the wavelength associated 
with the beat note Acj of the probe and control fields is 
comparable or shorter than L. For realistic values of L of 
a few centimeters, the term C2Q-D will have a significant 
impact if | Awj is of the order of a few GHz or larger. Note 
that polariton losses can be minimized by minimizing 
|Aa;|. This is obvious for £2^0 since it is proportional to 
Aw^. However, also the impact of £iPd increases with 
increasing values of | Aw | since cfcmax ~ | Aw | -I- cct^ . 

The total losses of dark-state polaritons can be calcu- 
lated from Eq. ([55[) . We find that the mean number of 



dark-state polaritons (Af) obeys 



inelastic polariton-polariton interactions. We obtain [ic 



dt{M)=-^ 



^2 



Auj' 



1 

W) 



k 



^fc(V'IV'fe)j(AA>, 
(64) 



fiQB = -iHcffgD + iQBHlff+Igu + fiCigB + fiC2gn, (66) 



where Af = '^j, V'lV'fe is the polariton number operator. 
The soUd Hne in Fig. |4] shows the losses of polaritons 
according to Eq. (|M|) where {iplfpk) was calculated via 
the numerical integration of Maxwell-Bloch equations, 
see Appendix |Dj On the other hand, the number of dark- 
state polaritons is proportional to the electromagnetic 
field intensity, 



(AA)(x dz{\g+f + \g^\') 



(65) 



where Q^ (C/_ ) is the Rabi frequency of the classical probe 
field propagating in the positive (negative) z direction. In 
order to test Eq. (|M)) . we evaluate the right-hand side of 
Eq. (|65p as a function of time from a numerical integra- 
tion of Maxwell-Bloch equations. The result is shown as 
the dotted line in Fig. |4] and in excellent agreement with 
the findings of Eq. (|M| . Finally, the dashed line in Fig.|3] 
represents the polariton losses if the term £2 Pd were ne- 
glected and shows that £2^0 contributes significantly. 
Note that the parameters of the first stationary light ex- 
periment [23 indicate a ratio of |Aaj|/(cfc„iax) ~ 0.75, 
which is even larger than the value of | Aaj/(cfcmax) ~ 0.25 
chosen in Fig. HI 

We emphasize that the loss term £2 go arises only in 
the presence of two counter-propagating control fields. 
In this case, the photonic component Ak in Eq. p^ 
of the dark-state polaritons tpk is comprised of counter- 
propagating probe field modes that are grouped around 
the wave numbers ztkc of the control field rather than 
the probe field, see Sec. IIIIBI It follows that the total 
Hamiltonian H — Ho+H\ + Hj^i^ in Sec.|n]does not pos- 
sess a true dark state for Aw 7^ 0. Even the dark-state 
polaritons "0^ in the fc = mode experience a coupling to 
bright-state polaritons (f>k and thus decay. This mecha- 
nism is at the heart of Eq. (j62p that describes the loss of 
dark-state polaritons in the fc = mode. The situation is 
different in a standard FIT configuration where only one 
pair of co-propagating probe and control fields is present, 
and the loss of dark-state polaritons is described by CigD 
only. 



B. STATIONARY LIGHT WITH 
TWO-PARTICLE INTERACTION 



Next we restate the full master equation (|^5|) in terms 
of the field operators tp{z) defined in Eq. ([M)) . In addi- 
tion to the terms discussed in Sec. IIVAI we have to take 
into account all contributions proportional to the cou- 
pling constant (72 in Eq. (HSI) that account for elastic and 



where i^eff is a non-hermitian Hamiltonian, 



ffcff - ^3 + I / dzV'tV' , 

^ Jo 



(67) 



and H3, Cigu and £2^0 are defined in Eqs. ([57]) . ((5^ 
and (|60p . respectively. The parameter 



2hLgj cos^ e 



^ A - cos2 6* Aw -I- 2742/2 
is the complex coupling constant, and 

L 

IgB = -Im(,g) / dzip"^ g-otp^"^ ■ 



(68) 



(69) 



The term proportional to g in Eq. (j67p and Xgu in 
Eq. (|69)) account for elastic and inelastic two-particle in- 
teractions that originate from the coupling of dark-state 
polaritons to the excited state |4). More precisely, the 
real part of g gives rise to a hermitian contribution to 
Hcs that accounts for elastic two-particle collisions. On 
the other hand, the imaginary part of g together with 
TgD gives rise to a two-particle loss term that can be 
written in Lindblad form as lm.{g/2)'D[ip'^]. 

The master equation ((66|) is equivalent to Eq. (j45|) and 
describes a one-dimensional system of bosons with effec- 
tive mass TTT-eff that experience elastic and inelastic two- 
particle interactions. Except for the two loss terms £iPd 
and £2^0, Eq. (|66)) can be identified with the dissipative 
Lieb-Liniger model discussed in the next Section. 



C. DISSIPATIVE LIEB-LINIGER MODEL 

The original Lieb-Liniger model [53 established in 
1963 describes bosons in one dimension that experience 
a repulsive contact interaction. In the limit of strong 
interactions, the bosons can enter the regime of a Tonks- 
Girardeau gas [5J| where they behave with respect to 
many observables as if they were fermions. Recently, it 
was shown [i^l that the original Lieb-Liniger model can 
be generalized to systems where the bosons experience a 
contact interaction with complex coupling constant, i.e., 
they undergo elastic or inelastic two- par ticle interactions. 
This dissipative Lieb-Liniger model [i^ shows that even 
a purely dissipative interaction effectively results in a re- 
pulsion and produces a Tonks-Girardeau gas in the limit 
of strong interactions. 

The master equation (|66p can be identified with the 
dissipative Lieb-Liniger model provided that the loss 
terms £i£>d and £2^0 are negligible. In the following 
we specify the conditions that justify this approxima- 
tion and assume that Ao;^ is small enough such that the 
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impact of C2QD is small as compared to CiQy), see Sec- 
tion llV Al On the other hand, the diffusion term CiQy^ is 
negligible if two conditions are met. First, the dynamics 
induced by the kinetic energy term proportional to mcs 
in Eq. (|67p must be fast as compared to the inverse de- 
cay rate of polaritons introduced by CiQd. This can be 
achieved if we set \5\ ^ F. Second, losses due to £igD 
must be negligible which imposes a limit on the maximal 
evolution time 



tn 



< 



1 1 



2nl 



Fcos2 6'c2fc2^^^' 



(70) 



Note that t,nax ^ 1/F is much larger than the lifetime of 
the excited state |3). Under these conditions, the master 
equation ((66)) reduces to 



hg-D 



-iHcftgu + igoHlg + Igu 



(71) 



and can be identified with the generalized Lieb-Liniger 
model [49[ for a one-dimensional system of bosons with 
mass rrieff and complex interaction parameter g. All fea- 
tures of the Lieb-Liniger model [4^, [SJ] are characterized 
by a single, dimensionless parameter 



G 



ft27Vph/£' 



(72) 



where TVph is the number of photons in the pulse. In the 
strongly correlated regime \G\ ^ 1, the interaction be- 
tween the particles creates a Tonks- Girardeau (TG) gas 
where polaritons behave like impenetrable hard-core par- 
ticles that never occupy the same position. Depending 
on the sign of the detuning 6 and A, the elastic inter- 
action between the polaritons can be either attractive or 
repulsive. The dissipative component of the interaction 
is negligible for A ^ 742. In this case, the preparation of 
a TG gas of polaritons with repulsion can be achieved if 
(5 A < [i2l- Note that the interaction becomes attrac- 
tive if i5A > which opens up the possibility to enter 
the super Tonks- Girardeau regime for polaritons [55l457| . 
Since the coupling constant g in Eq. (j68p is maximal for 
A = 0, the Lieb-Liniger parameter \G\ and hence the 
induced correlations are maximal for purely dissipative 
interactions 1 131. 



D. OTHER REALIZATIONS 

The master equation for dark-state polaritons in 
Sec. II VI was derived under the assumption that the level 
scheme of each atom is given by Fig. [21 Here we point 
out that stationary light and interacting dark-state po- 
laritons can be realized as well with the level scheme 
in Fig. [5] that was suggested in [ll|, [231 ■ Moreover, a 
straightforward modification of the formalism described 
in Sec. IIIII demonstrates that the level scheme of Fig. \5\ 
leads to the same master equation (|45|) for dark-state po- 
laritons as the configuration in Fig. [51 On the contrary, 



3) 



731 + 




* |1) 



FIG. 5. (Color online) Alternative level scheme that gives rise 
to the master equation for polaritons described in Sec. lIVI All 
symbols are defined in the same way as in Fig. [2] 



each system displays characteristic advantages and dis- 
advantages that we discuss now. 

The major difference between the two configurations 
is that the level scheme in Fig. [5l creates stationary light 
via the double- A system formed by states |1), |2), |3) 
and 1 3). Here each probe field interacts only with the 
co-propagating control field, and thus no fast oscillating 
spin coherences Xk {m) for m y^ are produced |58| . This 
feature is a significant advantage of the system shown 
in Fig. [5l since it implies that the master equation (|45)) 
remains valid for ultracold atoms or stationary atoms 
where the condition (1531) cannot be fulfilled. 



On the other hand, the implementation of the configu- 
ration in Fig. [SI comes along with difficulties that do not 
occur in the case of the level scheme in Fig. [2l First, we 
note that the transitions in the A configuration of Fig. [21 
can be selected by polarization, even if the ground states 
|1) and 1 2) are degenerate. This is not the case for the 
level scheme in Fig. [51 where the level splitting between 
the ground states |1) and |2) must be large enough such 
that the transitions |2) o |3), |1) o |3) and |2) o |3), 
|1) <-^ 1 3) can be addressed independently. Consequently, 
the parameter \Auj\ must be significantly larger than the 
Rabi frequencies r2± of the control fields which leads to 
additional losses, see Sec. [iVl Second, we point out that 
the initial polariton pulse in k space is centered around 
Auj/c if it is prepared via the slowing and stopping of a 
probe pulse. The reason for this is that the wave num- 
bers k of stationary-light polaritons have to be grouped 
around the wavenumbcrs ±fcc of the control fields, see 
Sec. nil Bl If the kinetic energy of the stationary-light 
polaritons is different from zero (i.e., S j^ 0) and in the 
case of the level scheme of Fig. [5l the mandatory choice 
of Acj y^ will lead to a moving polariton pulse even if 
the control fields have the same intensity. 
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V. SUMMARY 

In this paper wc introduced a technique for the de- 
scription of Ught-matter interactions under conditions of 
EIT. More specifically, we described a general method 
for the derivation of a master equation for dark-state po- 
laritons. In contrast to the standard description [la, [T^l 
based on a Heisenberg-Langevin approach, our master 
equation facilitates the treatment of polariton losses. 
This achievement allows us to model general polariton- 
polariton interactions that may be conservative, dissipa- 
tive or a mixture of both. In particular, the master equa- 
tion approach enables us to study dissipation-induced 
correlations j4^, [50| that are promising in the quest for 
highly correlated systems since they can be considerably 
stronger than their conservative counterparts [Tjl. For 
the illustration of our technique we use the example of 
stationary light polaritons that experience a conservative 
or dissipative interaction. The resulting master equation 
is discussed in various limiting cases. For stationary-light 
polaritons we find an additional loss mechanism that was 
overlooked so far. It is related to the fact that the to- 
tal Hamiltonian of the system does not possess a true 
dark state if the ground states are non-degenerate. In 
particular, polariton losses in level configurations with 
non-degenerate ground states can be a multiple of those 
in level schemes with degenerate ground states. Further- 
more, we specified conditions that allow us to reduce the 
full master equation for dark-state polaritons to the dis- 
sipative Lieb-Liniger model [40|. Finally, we discussed 
the atomic level scheme in Fig. [S]that leads to the same 
master equation for dark-state polaritons as the one in 
Fig. [2] and compared the advantages and disadvantages 
of both configurations. 

At the heart of our approach is a mapping of the full 
system dynamics to a conceptually simple system of cou- 
pled bosonic modes. This mapping could be the starting 
point for future studies of EIT systems beyond a Marko- 
vian master equation for dark-state polaritons. 



and that coincides approximately with X in the subspace 
"Hfe- The latter condition implies that Xq and X must 
necessarily obey the same commutation relations with 



the creation operators 0] {Oi G O) in 'Hfe, 



X,0 



■Hf 



Xo,0 



■Hf 



(Al) 



Furthermore, we have Xo\s} = X\s) for an arbitrary 
state \s) G 'Hfe- In all situations considered below, the 
operator X annihilates the vacuum state, XjO) = 0. In 
addition to Eq. (|A1[) . we thus require 



Xo\0)^0. 



(A2) 



The two conditions in Eqs. (jAip and (|A2p are sufficient 
to determine Xq since they guarantee that the matrix 
elements oi X and Xq are identical in the subspace "Hfe. 
In order to see this, we evaluate the action of X applied 
to the states \{ni, . . . , hm}) that span Hfe, 



X\{ni,. ..,nM}) = 



M 



^^Il-hiol) 



|0) . (A3) 



Here we employed the definition ([TT|) and the relation 
X\0) = 0. The recursive application of the identity 
[A, BC] = [A, B]C + B[A, C] allows us to write the com- 
mutator on the right-hand side of Eq. (jA3p as a sum 
of terms where only the commutator of X with one of 
the creation operators O] appears. Therefore, Eqs. (|A1[) 
and (|A2[) guarantee that the matrix elements of X and 
Xq are identical in the subspace Hfe- 

As an example, we discuss the representation of 



^N 



J2u=i ^22 in terms of the operators Oi € O. The only 



^M=i 



non-vanishing commutators are 
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Appendix A: REPRESENTATION OF 
OPERATORS IN Hfe 

Here we show how the master equation Eq. ([2|) can be 
expressed in terms of bosonic creation and annihilation 
operators if the system dynamics is restricted to the sub- 
space Hfe- To this end, let X be an operator acting on 
the total state space Htot- In the following, we describe 
a procedure that allows one to construct an operator Xq 
that comprises of bosonic operators Oi € O [sec Eq. 



iip) 



■ JV 
/ . ^22 '^'■k 



Xlim) 



Xlim), 



(A4) 



and we have X]u=i ^22 |0) 



N aMi 



0. According to 



Eqs. (|A1|) . (jA2p and with the bosonic commutation rela- 
tions obeyed by the operators Xk{m), the representation 

of E!Li ^22^ in ^FE is given by 



JV 

E^22^ - 5]5]4(m)X,(m). (A5) 



ij.=i 



m k 



The representation of the remaining operators that ap- 
pear in the Hamiltonian H in Eq. ^ can be found in a 
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similar way, the result is 

N 



fi—l -rn k 

N 

E4^' = EE^l(-)/fcM' 

N 



TV 

E^: 



m fe 



32 '^ 



^^X,(m)ift(,„_l), 



(A6) 

(A7) 
(A8) 
(A9) 



m k 



N 



TV 

E^^ 



it^)„i{k^+k)z^ 



42 



^^Xp_fe(m)/t(m+l), (AlO) 



?7l p 



M„-i(k^-k)z^ 
42 6 



^^X,_,(™)/t(m-l). 



(All) 



These relations together with the inverse relations of 
Eqs. (fT4)) and (|T5|) allow us to find the representation 
of Hq, Ha and i^NL in 'Hfe, 



(A12) 
(A13) 



m A: 



Ho = - fiAc^y^ (^I^fc + -Di-DA;) + 2^1 sin (^ cos (^ ^ Wfc (^[.^'fc + A-^^I 

- fe ^ ^ X^ (m)Xfe (m) - M ^ ^ if t (^)^^ (^) _ ;,(^ + g) ^ ^ jt (^) j^ (, 

rn k m k 

- fi(sin2 (^ - cos^ (^) Y^ LUk {olDk - AIAj^ , 

k 

Ha = -fif7osin6'^|Afe [i7^(l) sin (^ + FJ;(-1) cost 

k 

— hflo COS 2_, Ey "''^'^ ("^) "^fe (m + 1 ) sin (/? + i/| (m— 1 ) cos i^ 

m k 

Hnl = - ^92 E '^'^ E E ^p-kim) [ilim + 1) sin ip + ll{m - 1) cos ^] 

k m p 

-%2E-^'=EE^P-'=("^) [/^(TO + l)cosv3-/^(m- l)sini^] +h.c.. 



up 



+ Dk 



H, 



t(l)cos^-i7^(-l)sin^]} (A14) 



h.c. 



(A15) 



m p 



For the full representation of the master equation in "Hfe it remains to transform the decay term C^qd, see Eq. ([5]) 
The representation of the terms that describe the decay of the excited states in Eq. ([S]) can be found via Eqs. (jA6l 
and (|A7|) . Since the super-operator C^g must preserve the trace of the density operator, we find 

-^7^ = - ^ E E [Hl{ni)Hk{mye + gHl{m)Hkim) - 2Hk{m)gHl{m) 



m k 



(A16) 



742 

2 



E E [ll{'nn)Ik{m)~g + gll{m)Ik{m) - 2Ik{'m)gll{m) 



m k 



where F — 731 + 732 is the full decay rate of the excited state |3). 
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Appendix B: DEFINITIONS 

The master equation ([55]) is obtained from ([25)1 for 0+ = fi_ and if we exchange the set of operators 
{Ak,Xk{0),Hkil),Hki-l)Jkil)Jki-l)} by the new set {^k,^k,Pk,Qk,Uk,Vk} according to 



V-fc = Afecos6'-Xfe(0)sin6', 
Pk^[Hk{l) + Hki-l)]/V2, 
f/fc = [/fc(l) + /fc(-l)]/\/2, 



(jjk = Afe sin 61 + Xfe(O) cos 6*, 

Qk^[Hkil)-Hk{-l)]/V2, 
Vk = [Ikil)-Iki-l)]/V2. 



(Bl) 
(B2) 
(B3) 



Since the old and new operators are related by a unitary transformation, the new operators {"tpk, (j^k, Pk, Qk, Uk, Vk} 
obey bosonic commutation relations. If the inverse relations of Eqs. ()Bip - (jB3|) are plugged into Eqs. (jA12[) . ()A14|) 
(|A15|) and (|A16p . we obtain the master equation ([35]). Furthermore, we employ the condition in Eq. (j37p and the 
definition Ag = A — cot^ OAuj. The bath Hamiltonian in Eq. (P5|) is comprised of four parts, 



H^=Hi"^+Hi'^+Hi'^+H. 



r(3) 
B ' 



where 



(B4) 



<) ^-hAu;J2 [DiDk + 4<Pk] - hS Y, [PlPk + QlQk) - hAg ^ [ulUk + V^Vk 



hY^oJk (4>lDk + Dlcbk) - hno J2 {pUk + Pk4>l + QlOk + olQk 

k k 

fe ^ ^ xl{m)Xkim) - M ^ ^ Hl{m)Hk{m) -KAeY^Yl lli^)h{-m), 



(B5) 



ff, 



(1) 



-Mir, con 6 
2 



^ [Xk{-2) {Pl Ql) + Xk{2) [Pl + Ql 



^fif2ocos0^[ Y^ Xk{m)Hl{m + l)+ ^ Xk{m)Hl{m-l) 



V2 



+ h.c., 



(B6) 



F^2) = -hg^ cos eY \ul 



Y <Pp<Pk^P + ^k Y P'p^k-p 



h.c. 



(B7) 



4'^ = - ^f^m E E [{^l - ^k) Xk-pi~2) {^, + D,) + (Ul + V^) Xk-pi2) (0p - D,) 



k p 



4^52EE[ E li{m + mk^pim){^p + Dp)+ E ll{m-l)Xk^pi^)i(t>p-Dp 



(B8) 



+ h.c.. 



The dominant contribution to the bath Hamiltonian Hb is represented by Hj^ . The term H^ accounts for the 

(2) (3) 

modification of stationary light due to the fast oscillating spin coherences. Furthermore, Hj^ and i/g arise from the 
coupling of bath excitations to the transition |4) -r-> |2). We write the interaction Hamiltonian in Eq. ([55]) as a sum 
of three parts 



V = y(°) + ^(1) + T/(2) 



where 



(B9) 



1^(°) = - Hcos9AuY {d'Pk + ^k(f>i) +ncoseY^k [i^lDk + Dlij^ 

k k 

+ hg2 cos E ( ^l E ^P^'^-P + ^k E ^P^i-p ) ' 



(BIO) 



y(l) 



=fi52 J2 {^l E <^W'fe-P + Vk E ^p^'^-p) + i^-c- 
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^(2) 



lhg2CoseY,T.[{ul - ^/l) ^P^^-p(-2) + {Ul + V^) ^PpXk-p{2) 



k p 



-=%2COs6'^^| ^ 4(m + l)V'pXfc_p(TO)+ ^ ll{m - l)ippXk-pir 



^/2 



fe p 



+ h.c. 



(B12) 
(B13) 



The term V^^' describes the couphng of one or two dark-state polaritons to one bath excitation. On the other hand, 
V^^'^ arises from the couphng of one dark-state polariton and one bath excitation to the transition |4) <-> |2). The 
remaining part T^^^^ accounts for the couphng of a dark-state polariton and a fast spin coherence to the transition 



B);; 



1 4) o 1 2). Finally, the decay term C\ ' g in Eq 



reads 



4""'^ = - ^ E {pIP''~^ + bpIp^ - 2Pk~QPl) - ^ E (qIQ^^ + ^Ql^fc - 2gfeeQt 

fc fe 

' ^ E E [Hl{m)Hk{m.)g+ gHl{m)Hkim) - 2Hkim)gHl{m) 

- ? E {'^l^'^^ + -gulUk - 2Ukgul) -^Y. {Vkyke + gV^v^ - 2V,eV^) 
fe fe 

-^ E J2[4i^''^^k{m)g + gll{m)h{m)^2him)gll{rr 



(B14) 



Appendix C: BATH CORRELATION FUNCTIONS 



Here we outline the calculation of the bath correlation 
functions 



TrB{B,e^^"S]eB} 



(CI) 



that enter the master equation for dark-state polaritons 
via Eqs. (^0)) and (01]). In the following, we approximate 
the Liouvillian Cb by 



C 



(0) 



gB = -j:[H, 



(0) 
B : 



9d\ 



4«) 



QD, 



(C2) 



which amounts to neglect H^ , H^ and H^ in the ex- 
pression for Hb- In a first step, we show how the bath 
dynamics can be solved exactly with respect to the super- 
operator in Eq. (jC2p , and then we specify the conditions 



(1) 



r(2) 



r(3) 



that allow us to neglect H^ , H^ and iJg . The sim 



plified bath dynamics according to Eq. (|C2p justifies to 
replace V in the commutator [V, . . .] of Eq. (|4ip by V^^\ 
see Sec. IIIIDI Furthermore, we argue that the domi- 
nant contribution to Eq. (|44| stems from the interaction 
Hamiltonian V^^\ while F'^' and F'^' can be neglected. 
At the end of this section, we show that F'^' is indeed 
negligible if the dynamics of the system is restricted to 
the subspace "Hfe- In addition, all correlation functions 
in Eq. (|Cip where either Bi, Bj or both stem from F'-^-' 
vanish. V^'^' may thus only contribute to higher-order 



I 

terms beyond the Born approximation, but this effect 
will be small in the slow light limit since V^"^' is propor- 
tional to cos 61. 

Here we only take into account the bath operators 
Bi £ {Dk,4>k, Uk} that appear in the interaction Hamil- 
tonian Vq in Eq. (|B10[) . In principle, all combinations 
of these bath operators can enter the correlation func- 
tions in Eq. (|44)) . Their evaluation can be accomplished 
if the correlation functions in Eq. (|C1[) are regarded as 
mean values of an operator Bi with respect to the time- 
dependent, non-hermitian operator X 



e^^-B]gB. 



\B,\\ =Tr, 



{b.x] . 



(C3) 



It follows that the equations of motion for these mean 
values are given by 



dt\\B,\\ = -^TrB ([i3„<']x) +TrB [b.C^^^X 

(C4) 
If we apply this result to the operator Uk , we find that the 
time evolution of ||[/fc|| is determined by a single equation 



^,\\Uk\\ = {^Ao-^)m\\ 



(C5) 



On the other hand, the mean values of \\Dk\\ and \\4>k\\ 
are coupled to ||Pfe|| and ||Qfe|| via the following set of 
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linear equations, 

dtUkW = tAuj\\cj)k\\ - iujkWDkW + tnoWPkW , 
dt\\Pk\\ = ii5-T/2)\\Pk\\+inoUkh 

dt\\D4 - iAu\\Dk\\ - iLUkWMl + *^o||Qfc|| , 
dtWQkW = {i5 - r/2) IIQfcll + inoWDkW . (C6) 

It follows that the only non- vanishing terms in Eq. pi)) 
are given by 

Sigu) =Aw' cos2 e^Ii (v-i^p^D - i'pQBi^l) (C7) 

k,p 

+ cos^ e'^l2UJkUjp (ipli^pgD - ^pgui^l] 

k,p 

+ gl cos" OY^I^Y, 

k,p k' ,p' 

X (v'tV'Lfc'V'p'V'p-p'fi'D - V'p'V'p-p'^DV'tV'Lfc') 

- cos^ 0Aw ^ a;fel4 yJjliJpQn ~ ippQui'l) 

k,p 

- cos^ 9Auj ^ Wpis f V'fcV'pfi'D -ippQ-Dipi) , 

k,p 

where the integrals over the bath correlation functions 
are defined as 



where M is a 4 x 4 matrix and 



Ji(fc,_p) = / drTiB \4'ke 



Cbti^ 



Qb\ 



v{k,p) = j 



oo 

Mk,p) = fdrTTB [Dke'^^^DlgB] 



CJO 

Mk,p) = fdrTTB [Uke^^^U^QB] 



t 
Ti{k,p) = J drTvB [Dke^^^(j)lgB] 



t 
Mk,p) = J drTTB [^ke^^^DlQB] 



We illustrate the evaluation of the integrals in Eqs. (jC8 
(|C12p using the example of 



(C8) 

(C9) 

(CIO) 

(Cll) 

(C12) 



Ii = / drUkW 



(C13) 



where the mean value \\(/)k\\ is taken with respect to X = 
(jr qb- The integral in the latter equation can be regarded 
as the Laplace transform of \\(l>k\\ evaluated at s = 0. In 
order to determine Ii , we write the system of differential 
equations (|C6|) in matrix form, 



y = iUkl\\Pk\\,\\Dkl\\Qk\\)- 



(C15) 



Since all mean values tend to zero for t — > oo due to the 
presence of the decay term C\ Qd, the Laplace transform 
y{s) of y{t) exists and Eq. (|C14p yields sy{s) — y{0) = 
My{s). In the limit s ^ 0, we thus obtain 



y(0) = -A/-iy(0), 



(C16) 



where y(0) represents y at time t — Q. Since the mean 
values are taken with respect to X = (j))pQB, we have 
y(0) = (1,0,0,0) and thus y(0) can be determined. Fi- 
nally, Xi can be identified with [y(0)]i, i.e., the first com- 
ponent of 1/(0) in Eq, 



(jC16p . The evaluation of the re- 
maining integrals follows the same route and yields 



Re[Zi]==Rc[Z2]«^<^(fc,p), 
Im[2:i] = Im[l2] ~ -—i^5{k,p) , 

Im[Z3]=^2+'|^/4'5(fc,p), 

I4 = Xg W . 



(C17) 

(C18) 

(C19) 

(C20) 
(C21) 



These simple expressions for the integrals represent an 
expansion of more complicated terms that holds if JIq is 
sufficiently large as compared to the detuning \5\ and the 
decay rates of the excited states [see Eq. ([5(7)) ]. In ad- 
dition. Aw and |wfe|niax rnust be at most of the order of 
r^o- If the expressions in Eqs. (jC17|) - (|C2ip are plugged 
into Eq. (jC7p . we obtain the final result for our mas- 
ter equation (1^5]) . The validity of the Markov approxi- 
mation requires that the decay of the bath functions in 
Eqs. (JC8P - (|C12|) is fast as compared to the change of the 
density operator introduced by these terms. Since the 
correlation functions decay on a timescale that is of the 
order of l/7y, the Markov approximation is justified if 
the conditions in Eqs. ([5T|) and (j52p are met. 

Next we specify the conditions that allow us to neglect 



r(i) 



r(2) 



r(3) 



r(i) 



dty = My , 



(C14) 



H\^', H^' and H^' . The Hamiltonian H^'' describes 
the modification of stationary light due to the fast os- 
cillating spin excitations Xk{m) {m 7^ 0). These exci- 
tations are washed out due to the motion of the atoms, 
and the corresponding decay rate Fpo depends on the 
temperature of the atomic cloud. Note that the decay of 
the slowly varying spin excitations Xk (0) is significantly 
smaller than Fpo since the relevant wavenumbers are sev- 
eral orders of magnitude smaller. We emphasize that 
a Markovian master equation for the dark-state polari- 
tons corresponding to the level scheme in Fig. [5] is only 
possible if Fpo is comparable to the decay rate of the 
excited states. More specifically, the Hamiltonian H^ 
alters the set of equations (jC6|) and introduces a cou- 
pling between \\Pk\\, \\Qk\\ and the fast spin coherences 
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|lXfc(±2)||. Due to the decay of ||Xfe(±2)||, the effective 
couphng between ||Pfc||, \\Qk\\ and ||Xfe(±2)|j is given by 
(rioCOs6')^/rFo if r^ocos^ ^ Fpo 153 ■ It foUows that 
the effect of ifg is ncghgible provided that the condi- 
tion in Eq. (|53p holds, which means that the effective 
couphng between ||Pfe||, \\Qk\\ and ||Xfe(±2)|| is neghgi- 
ble on the timescale l/F which represents the hfetime 
of excitations in modes {Dk, 4>k, Pk,Qk}- Note that the 
decay of fast oscihating coherences and condition ([55]) is 
not required in the case of the level scheme discussed in 
Sec. lrTDl 

The Hamiltonian iJg gives rise to a modification of 
Eq. (fC5|) that determines ||C/fc||, 



dt\m\ 



iAfl 



'f) 



\Uk\ 



ig2 cos 



'El 
p 



Ofe-pl 



(C22) 



If the probe field modes form a (quasi-)continuum, then 
the second term in Eq. (|C22p gives rise to an additional 
decay channel of excitations in the mode Uk- We find 
that the associated decay rate is at most given by 



Tu^TiD^ocos^9/T, 



(C23) 



where Fid = 92^/c is the decay rate of the excited state 
|4) into the fiber modes. It follows that the influence 

(2) 

of _ffg ' is negligible provided that F[/ is much smaller 
than 742, which can always be achieved in the slow-light 
regime. 

It remains to discuss the impact of H^ and V''^\ 
These terms are negligible since the physical processes 
described by them are off-resonant and therefore strongly 
suppressed. Formally, the latter result can be obtained 
via a rotating-wave type approximation if the operators 
(/)fc and Dk in H^ and V^^' are expressed in terms of new 
bosonic operators that diagonalize -ffg . The influence of 
these operators is found to be small if 52/^0 oc 1/A^ <C 1, 
where N is the total number of atoms. 

Finally, we note that we have verified the validity of the 
approximations discussed above by the numerical com- 
parison of our master equation with the results of the 
full dynamics for a single mode. 



Appendix D: MAXWELL-BLOCH EQUATIONS 

Here we outline the numerical integration of the cou- 
pled Maxwell-Bloch equations [60| for classical probe and 
control fields that interact with the A subsystem formed 
by states |1), |2) and |3). The density operator of a sin- 
gle atom at position z is denoted by _R, and the coher- 
ence i?3i = (3|i?|l) is written as the sum of two counter- 
propagating terms. 



?(+) 



?(-), 



R3i{z,t) ^ R^^>iz,t)e'''''' + Rl,'iz,t)e 






(Dl) 



We assume that the probe fields are weak such that we 
can set Rn w 1, and apply the secular approximation 
where we drop fast oscillating terms exp\±2ik„z]. The 
Bloch equations of the system are thus given by [2^ [61| 

ati?^f = irj,eT''^^i?2i + (iS - T/2)R^^ + iQ^ (D2) 



dtR2Y=ieR2Y+i^cW~''RV^ 



?(+) 






(D3) 



where 5+ {Q- ) is the Rabi frequency of the classical probe 
field that propagates in the positive (negative) z direc- 
tion, and fcp = Wp/c is the wave number corresponding 
to the central frequency Wp of the probe field. Note that 
Q± [z, t) are slowly varying functions of position and time. 
On the other hand, the Rabi frequency Q,c of the each 
control field is assumed to be position-independent but 
varies with time. The Bloch equations have to be solved 
consistently with Maxwell's equations that yield [2^, |61| 



^•^ ^F,\r -^ffi 7,(±) 



(D4) 



where we employed the slowly varying envelope approx- 
imation [601 • The set of equations (jD3[) and (jD4[) allows 
us to determine g± as well as the atomic variables R. 
Note that equivalent results without the secular approx- 
imation can be obtained |23| if the level scheme in Fig. |3] 
is employed. 

In the slow light limit, the expectation value of the 
polariton pulse in position space is directly proportional 
to the expectation value of the ground-state coherence. 



mz)) (X (i?i2>(z). 



(D5) 



It follows that (ipk) can be calculated from (i?i2)(z) via 
a Fourier transformation with respect to position, and 



we have (V'lV'fc) 

corresponds to a coherent state. 



KV'fc)!^ since a classical probe field 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, 2000). 

[2] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 



80, 885 (2008). 
[3] M. J. Hartmann, F. G. S. L. Brandao, and M. B. Plenio, 

Nat. Phys. 2, 849 (2006). 
[4] M. J. Hartmann, F. G. S. L. Brandao, and M. B. Plenio, 



17 



Laser & Photon. Rev. 2, 527 (2008). 
[5] D. Rossini and R. Fazio, Phys. Rev. Lett. 99, 186401 [36 

(2007). 
[6] D. G. Angelakis, M. F. Santos, and S. Bose, Phys. Rev. [37 

A 76, 031805(R) (2007). 
[7] D.Gerace, H. E. Tiireci, A. Imamoglu, V. Giovannetti, [38 

and R. Fazio, Nat. Phys. 5, 281 (2006). [39 

[8] M. Hafezi and D. E. Chang and V. Gritsev and E. Demler 

and M. Lukin, .arXiv:0907.5206. (2009) . [40 

[9] M. Hafezi and D. E. Chang and V. Gritsev and E. Demler 

and M. Lukin, ■arXiv:0911. 4766' (2009). [41 

[10] I. Carusotto, D. Gerace, H. E. Tureci, S. De Liberato, 

C. Ciuti, and A. Imamoglu, Phys. Rev. Lett. 103, 033601 [42 

(2009). 
[11] M. Fleischhauer, J. Otterbach, and R. G. Unanyan, Phys. 

Rev. Lett. 101, 163601 (2008). [43 

[12] D. E. Chang, V. Gritsev, G. Morigi, V. Vuletic, M. D. 

Lukin, and E. A. Demler, Nat. Phys. 4, 884 (2008). 
[13] M. Kiflner and M. J. Hartmann, Phys. Rev. A 81, [44 

021806(R) (2010). [45^ 

[14] M. J. Hartmann, F. G. S. L. Brandao, and M. B. Plenio, 

New. J. Phys. 10, 033011 (2008). 
[15] M. J. Hartmann, Phys. Rev. Lett. 104, 113601 (2010). [46 

[16] M. Fleischhauer and M. D. Lukin, Phys. Rev. Lett. 84, 

5094 (2000). 
[17] M. Fleischhauer and M. D. Lukin, Phys. Rev. A 65, [47 

022314 (2002). 
[18] M. D. Lukin, Rev. Mod. Phys. 75, 457 (2003). 
[19] M. Fleischhauer, A. Imamoglu, and J. P. Marangos, Rev. [ 

Mod. Phys. 77, 633 (2005). 
[20] L. V. Hau, S. E. Harris, Z. Dutton, and C. Behroozi, 

Nature 397, 594 (1999). [49 

[21] M. M. Kash, V. A. Sautenkov, A. S. Zibrov, L. HoUberg, 

G. R. Welch, M. D. Lukin, Y. Rostovtsev, E. S. Fry, and 

M. O. Scully, Phys. Rev. Lett. 82, 5229 (1999). [50 

[22] D. Budker, D. F. Kimball, S. M. Rochester, and V. V. 

Yashchuk, Phys. Rev. Lett. 83, 1767 (1999). 
[23] C. Liu, Z. Dutton, C. H. Behroozi, and L. V. Hau, Nature [51 

409, 490 (2001). 
[24] D. F. Phillips, A. Fleischhauer, A. Mair, R. L. 

Walsworth, and M. D. Lukin, Phys. Rev. Lett. 86, 783 [52 

(2001). 
[25] T. N. Dey and G. S. Agarwal, Phys. Rev. A 67, 033813 [53 

(2003). [54 

[26] A. Andre and M. D. Lukin, Phys. Rev. Lett. 89, 143602 [55 

(2002). 
[27] M. Bajcsy, A. S. Zibrov, and M. D. Lukin, Nature 426, [56 

638 (2003). 
[28] F. E. Zimmer, A. Andre, M. D. Lukin, and M. Fleis- [57 

chhauer. Opt. Commun. 264, 441 (2006). 
[29] F. E. Zimmer, J. Otterbach, R. G. Unanyan, B. W. Shore, 

and M. Fleischhauer, Phys. Rev. A 77, 063823 (2008). [58 

[30] S. A. Moiseev and B. S. Ham, Phys. Rev. A 71, 053802 

(2005). 
[31] S. A. Moiseev and B. S. Ham, Phys. Rev. A 73, 033812 [59 

(2006). 
[32] Y.-W. Lin, W.-T. Liao, T. Peters, H.-C. Chou, J.-S. 

Wang, H.-W. Cho, P.-C. Kuan, and I. A. Yu, Phys. Rev. [60 

Lett. 102, 213601 (2009). 
[33] G. Nikoghosyan and M. Fleischhauer, Phys. Rev. A 80, [61 

013818 (2009). 
[34] H. Schmidt and A. Imamoglu, Opt. Lett. 21, 1936 (1996). 
[35] A. Imamoglu, H. Schmidt, G. Woods, and M. Deutsch, 



Phys. Rev. Lett. 79, 1467 (1997). 

S. E. Harris and Y. Yamamoto, Phys. Rev. Lett. 81, 3611 

(1998). 

S. E. Harris and L. V. Hau, Phys. Rev. Lett. 82, 4611 

(1999). 

H. Kang and Y. Zhu, Phys. Rev. Lett. 91, 093601 (2003). 

D. A. Braje, V. Balic, S. Coda, G. Y. Yin, and S. E. 
Harris, Phys. Rev. Lett. 93, 183601 (2004). 

M. J. Hartmann and M. B. Plenio, Phys. Rev. Lett. 99, 

103601 (2007). 

A. Andre, M. Bajcsy, A. S. Zibrov, and M. D. Lukin, 

Phys. Rev. Lett. 94, 063902 (2005). 

T. Aoki, B. Dayan, E. Wilcut, W.P.Bowen, A. Parkins, 

H. Kimble, T. Kippenberg, and K. Vahala, Nature 443, 

671 (2006). 

A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. 

Huang, J. Majer, S. Kumar, S. M. Girvin, and R. J. 

Schoelkopf, Nature 431, 162 (2004). 

M. Leib and M. J. Hartmann, arXi v:1006.2935l ('20101. 
K. Hennessy, A. Badolato, M. Winger, D. Gerace, 
M. Atature, S. Guide, S. Fait, E. L. Hu, and 
A. Imamoglu, Nature 445, 896 (2007). 
M. Trupke, J. Goldwin, B. Darquie, G. Dutier, S. Eriks- 
son, J. Ashmore, and E. A. Hinds, Phys. Rev. Lett. 99, 
063601 (2007). 

M. Bajcsy, S. Hofferberth, V. Balic, T. Peyronel, 
M. Hafezi, A. S. Zibrov, V. Vuletic, and M. D. Lukin, 
Phys. Rev. Lett. 102, 203902 (2009). 

E. Vetsch and D. Reitz and G. Sague and R. Schmidt and 
S. T. Dawkins and A. Rauschenbeutel, larXiv:0912. 11791 
(2009). 

S. Diirr, J. J. Garcia-RipoU, N. Syassen, D. M. Bauer, 

M. Lettner, J. I. Cirac, and G. Rempe, Phys. Rev. A 79, 

023614 (2009). 

N. Syassen, D. M. Bauer, M. Lettner, T. Volz, D. Dietze, 

J. J. Garcia-RipoU, J. I. Cirac, G. Rempe, and S. Diirr, 

Science 320, 1329 (2008). 

H.-P. Breuer and F. Petruccione, The Theory of Open 

Quantum Systems (Oxford University Press, Oxford, 

2006). 

C. Liu, Z. Dutton, C. H. Behroozi, and L. V. Hau, Nature 

409, 490 (2001). 

E. H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963). 

M. Girardeau, J. Math. Phys. 1, 516 (1960). 

G. E. Astrakharchik, J. Boronat, J. CasuUeras, and 

S. Giorgini, Phys. Rev. Lett. 95, 190407 (2005). 

M. T. Batchelor, M. Bortz, X. W. Guan, and N. Oelkers, 

J. Stat. Mech. L, 10001 (2005). 

E. Haller, M. Gustavsson, M. J. Mark, J. G. Danzl, 

R. Hart, G. Pupillo, and H.-C. Nagerl, Science 325, 1224 

(2009). 

Note that fast oscillating spin coherences may still be 

produced via the |2) <->■ |4) transition, but this process is 

off-resonant and thus suppressed (see Appendix [Cj. 

C. Cohen- Tannoudji, J. Dupont-Roc, G. Grynberg, 

Atom-Photon Interactions (Wiley, New York, 1992), Sec. 

in.c.3. 

M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- 
bridge University Press, Cambridge, 1997). 
M. Kiffner and T. N. Dey, Phys. Rev. A 79, 023829 
(2009). 



